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ABSTRACT 


The determination of overpressure yield from internal 
blast has been the subject of prior investigations. Models 
have been developed to predict explcsion thermodynamic param- 
eters, based on equations of state, mass and energy conserva- 
tion, and rate laws, utilizing high speed digital computers. 
Typically, thermodynamic equilibrium was assumed, an equilib- 
rium temperature was determined, ideal gas law behavior was 
ascribed to product gases, and a pressure was computed. 
Reports of the results of these studies have typically con- 
tained a disclaimer that radiation effects have been ignored. 
In this investigation, a blackbody radiator model was developed 
and, utilizing the parameters reported in the literature for a 
number of explosives of common interest, a determination was 
made of the effect of thermal radiation as an attenuation 
mechanism on explosive energy and the existence of thermody- 
namic equilibrium. Energy losses of less than .5% were calcu- 
lated and thermodynamic equilibrium was demonstrated through 


statistical considerations. 
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I. INTRODUCTION 


Internal blast is the term applied to the concussion 
effect resulting from highly exothermic reactions occurring 
in a confined space. Examples of internal blast are dust 
explosions in grain elevators and coal mines; gas pocket 
explosions in buildings and underground systems such as 
sewers and mines; and the detonation of high explosives, 
either accidentiy or through proper fuse function, ina 
cherge hole in a mine ore face, or in a shipboard ordnance 
stores magazine or other compartment. Substantial effort 
has been devoted to the study of explcsive phenomena during 
the past century; internal blast investigation has been 
particularly intense during the past two decades. 

Various computer codes have been developed to predict 
the thermodynamic properties of the detonation products 
generated by explosive materials. Typically, these codes 
have been used to determine a product composition and condi- 
tion based on chemical reaction rate equations and other 
thermodynamic and conservation considerations applied to a 
system of equations of state, such as the STRETCH code, which 
is based on the Becker-Kistiakowsky-Wilson (B-K-W) equations 
of state, reported by Mader of the Los Alamos Scientific 
Laboratory [Ref. 1). 

Other work has been very specifically directed to the 
determination of internal blast overpressure yield. Previous 


investigations were conducted by G. F. Kinney, R. A. Reinhardt 





and K. J. Graham of the Naval Postgraduate School, working 

in conjunction with personnel from the Naval Weapons Center, 
China Lake, California, on determination of peak overpressure 
yield of various explesives and fuels, with and without the 
incorporation of blast-augmenting reactive metals [Ref. 2, 

3 and 4]. The technique for determining internal blast 
overpressure used by Kinney, Sewell and Graham and by 
Reinhardt is to solve a system of simultaneous equilibrium 
equations to determine a product composition, based on a 
selected trial temperature. The sum of the internal energies 
of each product gas species at this trial temperature is 

then computed, using a polynomial form for C., for each 
species. The value of the sum is compared to the previously 
determined value of the maximum energy available for release 
in the explosion. The value of the trial temperature is 
varied until energy balance is achieved. This temperature is 
an equilibrium flame temperature. Ideal gas law behavior is 
then ascribed to the mass of product gases, and a peak over- 
pressure is computed. 

A second major effort was undertaken at the Naval Ordnance 
Laboratory, White Oak, Maryland, with a focus on the produc- 
tion of a computer program capable of predicting damage 
mechanism due to internal blast [Ref. 5]. A common disclaimer 
made in the reports of these investigations has been that the 
effects of thermal radiation as an energy loss mechanism can 
be ignored. Additionally, an assumption of thermodynamic 


equilibrium was made in the determination of flame temperature. 





The validation of the assumption of thermodynamic equilibrium 
establishment in the detonation reaction zone and the magni- 
tude of the thermal radiation losses are the subject of this 
investigation. Examination of open literature reports of 
investigations conducted on luminosity of explosions was made 
to determine if thermal radiation energetics in conventional 
high explosive detonations had been quantified. 

Inconsistencies of values which various sources cite for 
the same parameters is noted in examining the literature 
reporting explosives experimental results. A remark by 
Rudlin [Ref. 6] to the effect that statistical deviations 
within explosives data sets ostensibly reporting replicate 
measurements average about 10 to 15%, and that agreement to 
within 5% was unusually good, provides a basis for the order 
of magnitude and two digit precision used in the model 
development for this investigation. The best experimental 
or computational values available are used, the best being 
defined as having been cited in more than one source, or 
being internally consistent with other data cited in more 
than one source. 

A model of a detonating explosive is developed to explcre 
the existence of thermodynamic equilibrium and to quantify 
expected radiant exitance from the detonation. This model 
is established as a blackbody, with size and temperature 
characteristics determined by the detonation parameters of 
the explosives considered, as found in the literature. Util- 


izing these parameters, a determination is made of the effects 








of thermal radiation as an attenuation mechanism on the explo- 
Sive energy available for generation of overpressure, and of 


the existence of thermodynamic equilibrium. 





II. DETONATION THEORY 


A brief discussion of the theory of detonation and terms 
used in detonation literature is considered appropriate prior 
to discussion of the specific model developed. Condensed 
explosives are solids which consist of materials or mixtures 
of materials which are capable of undergoing highly excthermic 
chemical reactions at extremely fast rates, producing gaseous 
and, usually, solid reaction prcducts at high pressures and 
temperatures. "CHNO" explosives are organic materials whose 
molecules contain atoms of carbon, hydrogen, nitrogen and 
oxygen. These complex molecules undergo a decomposition reac-~ 
tion, followed immediately by oxidation-reduction reactions 
which produce low molecular weight detonation products, such 
5 No and C. The exact quantity of each reaction 


product, per mole of explosive, is dependent on such factors 


as HO, 2cO, CO 


as the initial composition and density of the explosive, the 
fuel-to-air ratio and initial conditions of the environment. 
Detonation of an explosive results from rates of reaction 
which are so rapid that the heat evolved in the exothermic 
decomposition and subsequent oxidation-reduction processes 
cannot be carried away before there is an energy build-up of 
such magnitude that a shock front results. This shock front 
propagates through the explosive materiai; it is called the 
detonation wave, and propagates at a characteristic velocity 


of detonation, D. 


10 





The detonation shock front separates the unreacted explo- 
sive material from the material undergoing reaction. One 
description (the model adopted for this investigation) considers 
three areaS in an explosive undergoing detonation, separated by 
two planes, as depicted in Figure 1. In this description, the 
shock wave initiates the decomposition reaction of the complex 
"CHNO" molecule. Decomposition reactions are completed and 
oxidation-reduction reactions proceed to completicn in the 
reaction zone. When the "C-J plane" is reached, reaction 
products are at thermodynamic equilibrium, and are of composi- 
tion and condition predictable by appropriate equations of 
Sir arce . 

Part of the net energy in the reaction zone pumps the 
detonation wave to sustain it as it propagates through the 
explosive. The rest of the net reaction energy is taken up 
as internal energy of the reaction products. When the detona- 
tion wave reaches the surface of the explosive, it shocks the 
Surrounding medium, setting up the initial shock wave in the 
medium, and establishes a transient reflected shock wave 
back into the material, which results in a compression. Thus, 
the net energy released in the detonation is partitioned, with 
some gOing into the surrounding medium as the initial shock- 
wave, and the rest being retained in the explosive material. 

The values of specific heat of detonation, q, reported 
in the literature are cften based on constant volume bomb 
calorimetric measurements. Values of heat of detonation 


experimentally determined in calorimetric devices are lower 


id 
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than the maximum heat of detonation calculated for the enthalpy 
of the reactions transforming the explosive to its detonation 
products [Ref. 8]. The bomb calorimeter values are lower 

than calculated values because the products formed in the 

bomb differ from those in the Chapman-Jouguet (C-J) detonation 
state. The final products achieved in the bomb calorimeter 

are not the C-J products but, rather, reflect approximate 
equilibrium at low pressures. The relative amounts of H.,0O, 


2 


CO. and solid carbon formed in the reactions, as compared to 


2 


CO and H are determined by pressure dependent equilibrium 


9! 
equations. The formation of HO, CO. and solid carbon are 
preferred over the formation of CO and H. at the very high 
pressure of a C-J detonation state. The detonation products 
analyzed from bomb calorimeter tests, even when confined, 


show much higher concentrations of CO and H. than predicted, 


2 


at the expense of H.O, CO. and solid carbon [Ref. 8]. The 


2 Z 
enthalpy of the C-J detonation products is more negative than 
that for the bomb products. Since the object of this investi- 
gation is to determine the thermal radiant loss from the more 
energetic detonation state, calculated values are used [Ref. 
7), where applicable. 

A brief review of the pertinent points of the theories 
which have been advanced to explain explosive behavior pro- 
vides definition to the terminology and concepts presented 
above. Reference 8, chapter 6, contains a satisfactory devel- 


opment of the elementary theory of the plane wave detonation. 


Extension of the concepts presented there to a spherical wave 


ine 





is straightforward [Ref. 8, ch. 9]. The detonation wave, in 
the framework of elementary thermodynamic-hydrodynamic theory, 
is assumed to be a shock wave, treated as a mathematical dis- 
continuity, with full application of the general conservation 
laws which apply across any such discontinuity. As such, then, 
the detonation wave obeys the Rankine-Hugoniot equation for a 
shock wave: 


ed aa "(Dy i a) Ue 7 vow 


where the subscript o refers to the unreacted material and 

1 to the reaction products; e is the specific internal energy 
function, p is the pressure and v is the specific volume. 
This equation will be referred to as the Hugoniot. If the 
products are assumed to be in thermodynamic equilibrium, then 
their thermodynamic state is specified as 


ey) = e(p,/V,) 


where e(p,v) is some known function of the state variables. 
In terms of specific enthalpy, 
h =e + pv, 


the Hugoniot can be written in an alternative form 


Ay - Ae = “(Py = Po) (vy + Vo): 


A plot of the Hugoniot in the Py-Vv plane provides a frame- 


a 
work for definition of the Chapman-Jouguet (C-J) point (see 
Figure 2). 

Consider the point on the Hugoniot curve for which Vi 
1s equal to oe Corresponding to conversion of the explosive 


to its products at constant volume. The value of pressure 
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Py for this volume Vv on the Hugoniot, Pa is the minimum 
value of pressure for detonation, and defines the lower 
point (Do +V)) of the detonation limb of the Hugoniot. Mass 


and momentum conservation conditions yield the relationship 


which describes a straight line in the P)- Vv plane; it is 


Ht 
called a Rayleigh line, and is the chord connecting some 
point (Py Vy) to the point (DorV,) - The propagation velocity, 
D, 1S proportional to the square root of the negative slope 
of the Rayleigh line. One Rayleigh line passes through 
(DorV) and is tangent to the Hugoniot on the detonation limb. 
The point of tangency is called the Chapman-Jouguet point 
(Pa_z"fc_z) and the explosion at this point has several 
important properties. The C-J state detonation, from some 
specified initial state (Dov a). has a minimum value for the 
propagation velocity D, relative to the unreacted material, 
and a stationary minimum value of entropy. The Hugoniot 
curve is tangent to the isentrope through the C-J point, and 
therefore the Rayleigh line is also tangent to that same 
isentrope. The flow of the reaction products relative to 
the reaction front is sonic behind a C-J detonation. The 
Chapman-Jouguet hypothesis states that the detonation wave 
for which the wave flow proceeds at a sonic velocity is the 


one that actually occurs when the reaction wave is unsupported. 


subsequent arguments by Becker [Ref. 9] showed, through 
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entrophy considerations, that only the C-J detonation is 
stable. The C-J plane is the locus of C-J points, where the 
chemical reactions have come to thermodynamic equilibrium and 
the detonation properties predicted by the equations of state 


exist. 


Ny 





III. MODEL DEVELOPMENT 


A. GENERAL DESCRIPTION AND GEOMETRY 

This investigation is Ber ateataadlth the centrally initiated 
detonation of spheres of condensed organic explosives of 
"CHNO" composition both with and without aluminum. A standard 
mass of ten kilograms is modeled, contained ina sphere of 
about 10 to 15 centimeters radius, depending on the density 
of the explosive. The volume taken up by the initiator at the 
center is ignored as is any energy contribution from the 
initiator, since both are very small compared to those quanti- 
ties for the main charge. The explosive is detonated ina 
rigid container of irregular geometry, with an average distance 
from the charge to the wall of five meters, filled with normal 
air at standard conditions. When the detonation shock wave 
reaches the surface of the explosive, the initial shock wave 
is established in the surrounding air, there is a contraction 
of the surface of the explosive due to the reflected shock 
Bemeweeme tne reaction zone is the surface of the sphere. 
The shock wave displaces the air slightly away from the sur- 
face of the sphere and the compression of the surface itself 
results in a rarefaction; any chemical interaction of the 
medium with the material in the reaction zone is initially 
ignored, for reasons developed below. The reaction zone, at 
the surface, is referred to as the fireball in some contexts; 
after.the hot gases produced by the detonation reaction have 


begun expansion and encounter material from the surrounding 
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medium, the surface does indeed burn in some cases. This is 
called afterburn and He most evident in oxygen-poor explosives 
or those mixed with reactive metals. For purposes of this 
model, any medium-product reactions that occur do so after 

the beginning of the hot gas expansion and are considered to 
slow down the decay of the high temperature-high pressure 
conditions when volume expansion occurs. 

The surface of the explosive is modeled as a blackbody. 
This requires thermodynamic equilibrium, which is demonstrated 
below. Spherical geometry is selected for several reasons. 

An isotropic blackbody obeys the Stefan-Boltzmann law, anda 
spherical surface is isotropic. For most explosive configura- 
tions of interest, at some time after completion of detonation, 
the expanding mass of hot gases assumes an approximately 
spherical geometry. The true geometry of a real system of 
explosive and confining chamber can be transformed into a 
spherical coordination description. Spherical geometry 
greatly facilitates computations of radiant energies. Addi- 
tionally, the experimental results reported in the literature 
for many explosives of interest have been for spherical masses 
of explosive. 

The initiation of detonation of a spherical mass of con- 
densed high explosive from the center produces an expanding 
spherical detonation shock wave which propagates through the 
explosive. Although several sources have indicated a problem 
with maintaining a spherically diverging finite reaction zone 


of constant width [Ref. 8], because it would of necessity 


1 





decrease in width as it expands, Cook [Ref. 10] employs a 
spheriform shock front model, and photographic evidence for 
such an expansion was reported by Rudlin [Ref. 11]. A point 
of consideration to resolve this quandry is that the value 

of the propagation velocity, D, builds up to the stable C-J 
value more slowly in the spherical case than in the case of 
plane wave propagation through an end-initiated cylindrical 
charge, and is a function of the radial distance out from the 
center [Ref. 8, ch. 9]. For purposes of the model used in 
this investigation, the C-J conditions are considered to exist 
at the surface of the sphere at the instant of departure of 
the initial shock wave into the medium. A minimum size exists 
below which this model does not apply, due to not reaching C-J 
conditions by the time the detonation wave has reached the 
surface, and due to an unfavorably small mass-to-surface area 
ratio (since heat of detonation and, therefore, the energy 
evolved, is a function of the mass of exploSive and the radi- 
ant energy emitted is a function of the surface area and 


temperature). 


B. THERMODYNAMIC EQUILIBRIUM 

Existence of thermodynamic equilibrium in the detonation 
reaction zone during a sensible time period is an essential 
element of the foundation of assumptions upon which internal 
blast studies have been built. Sensible time is chosen to be 
that time required for the reaction zone to move through the 
detonating explosive a distance equal to the width of the 


reaction zone. Although an exact value for the width of a 


20 





reaction zone is dependent on a number of variables, such as 
initial conditions of the explosive (density, packing, configu- 
ration, temperature, etc.) and the chemical composition and 
experimental measurement is, at best, very difficult, values 
on the order of 0.1 millimeter have been reported, based on 
pressure distribution in the reaction zone determined by 
shockwave measurements in thin layers of inert material placed 
in contact with the end surface of the charge [Ref. 12]. In 
Ref. 10, Cook makes a correlation between the reaction zone 
width and an area of ionization which exists immediately be- 
hind the detonation front. Using this correlation, reaction 
zone widths on the order of centimeters are reported. The 
reaction zone does propagate through the explosive immediately 
behind the detonation wave, with essentially the same velocity, 
D. As a typical case, TNT is considered. (Appendix A is a 
glossary of explosive acronyms; the common names of the explo- 
Sives considered in this investigation are defined, and their 
compositions are given.) Reported reaction zone widths of 

2.5 cm [Ref. 10] and 0.1mm [Ref. 12] are used to give an 


6 


upper limit value for the sensible time of 3.7 x 10 -s and 


wer lanit yaluesof 1.5 x 107° 


s, based on a reported 
detonation velocity value of 6830 m/s [Ref. 1]. If a large 
number of collisions occur between molecules of the product 
gases during this sensible time, it is reasonable to conclude 
that thermodynamic equilibrium exists. 


Geometric collisional diameters of product species have 


been reported in the range of one to five Angstroms (1 to 


2. 








5 omen determined by viscosity, heat conductivity and 


diffusion experiments [Ref. 13] or molecular beam apparatus 
experiments [Ref. 14]. Values of such collision diameters 


for diatomic and polyatomic molecular species are found to 


Leer H, and 4.82 x tO meter 


CO., with many values reported in the range of 3.5 to 5.0 x 


ae for molecules of the species found in detonation 


be, for example, 2.47 x 10 


product gases. Order of magnitude agreement was found through- 
out the literature examined. 

The number of particles, their collisional cross sections 
and thermal velocities can be combined to give a measure of 
the average time between collisions, through the relationship 


= 1 
6 * Tov 


where n is defined as the number density of particles in the 
gas, v as the average velocity and o as the collisional 
cross section. Reference 15 contains a form of expression 
for v (the Maxwell mean speed); using the variables of this 
analysis, v = (8kT/mm) 7. Reference 16 cites an equivalent 
expression, and provides a method to obtain an estimate of 
cross section for collisions between molecules of different 
species if molecular radii are known: o = 7 (ro is rp)? 
where re is the radius of molecule species "a" and rh is 
that of molecule species "b." Variations of molecular diam- 
eters, and the large number of possible varieties of colli- 


sion events that could occur between the dozen species usu- 


ally found in product gases, call for a simplifying 
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approximation. An average molecular collision diameter of 


10 


4x 10m is selected, which results in a value for o of 


-~20 2 | See 
m roy (tin lar 


Borx 10 (ce 10” 
Composition of product gas for a large number of organic 
chemical high explosives is reported by Mader [Ref. 1]. An 
average molecular weight for these gases was determined to 
be 30 (28.8 for TNT of density 1.60 gm/cm*, as an example). 
The particle density, n, is determined by multiplying the 
Loschmidt number (2.69 x Nowe particles/cm’) by the detona- 
tion pressure in atmospheres and by the fraction 273/ 
(detonation temperature). For condensed TNT of density 
1.60 gm/cm>, Ref. 1 reports detonation pressure of 194 kbars 
and detonation temperature of 2971 K. Using these values 
and an average molecular weight of 28.2, n was calculated to 
be 4.73 x woos particles/cm’, v was calculated to be 1.49 
x 1 cm/s, and a value of tO was calculated to be 2.84 x 


tome. Comparison of te to the shortest sensible time 


Veo ret.5 x 10°° 


s demonstrates that a very large number 
of collisional events, on the order of millions, occur during 
the sensible time. Thermodynamic equilibrium thus exists in 


the reaction zone at the C-J plane. 


C. BLACKBODY APPROXIMATION 

The maximum radiant energy emitted by any isotropic body 
at thermodynamic equilibrium is, by definition, that radiated 
by a blackbody. A determination of the spectral energy given 
off by a blackbody at the detonation temperature for some 


explosive provides a measure of the maximum energy which 
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could be lost from the product gases through thermal radia- 
tion. The surface of a sphere of detonated explosive, immedi- 
ately after the departure of the shock wave into the 
surrounding medium, is modeled as a blackbody made up of hot 
gas which emits the spectrum characteristic of the tempera- 
ture of detonation (the C-J temperature). Justification of 
this model follows order of magnitude comparisons between the 
excited state lifetimes and the time between collisions in 

the gas. 

The lifetime of an excited state is defined as the aver- 
age time an atom spends in this state before making a transi- 
tion. Radiation can be given off by the atom when transition 
to a lower state occurs. Lifetimes for many states are on 
the order of nanoseconds lon s) [Ref. 17]. Comparison of 
typical excited state (radiator) lifetimes orom S) to the 
average time between collisions for the gas molecules in a 


15 


detonation (10 ~~s), as computed above, leads to the conclu- 


sion that the individual spectral line widths for the dozen 
or So species present in the detonation reaction zone will 

be collisionally broadened to the point of showing a continu- 
um. Consider, as an example, the case of the 10.2 micron 


CO, line. It has a center frequency of 2.94 x 101 7Hz. e@) 


2 2 
has a collision diameter of 4.84 x io) ae At a C-J tempera- 
ture of 3000 K and pressure of 2 x 10> atmospheres, a time 


Ibs. 


of collision of about 10 ~"°s is calculated. The line width 


(frequency spread) due to collisional broadening is about 
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‘A aie 1 
- To [Ref. 18]. 


In this example, a line width of about 3x Ge He results, as 
compared to a room temperature width of about 50 MHz. This 
collisional broadening smears the 10.2 micron line from 

about 10°Hz through 107° Hz. The line spectrum of an explosion 
is very rich [Ref. 19]; such collisional broadening is opera- 
ting on all species line widths to some comparable degree, 
Causing them to smear together into an apparent continuum. 
Within the excited state lifetime, hundreds of thousands of 
collisions occur. Thus, thermodynamic equilibrium and continu- 
um radiation are characteristic of the reaction zone, and 
application of a blackbody model is quite justified. 


The radiant exitance of a blackbody is given by the 


Stefan-Boltzmann law, 


2 ae 


4 cal cm “~ K ie 


mei where co = 1.36x10°"- 


Using this relationship, the energy emitted by the model black- 
body can be found if the surface area, temperature and time 
of radiation are known. 

The radiant surface of the detonation is the surface of a 
sphere occupied by a ten kilogram mass of condensed explosive, 
and a function of the density of the explosive. Different 
loading techniques produce different densities for the same 
explosive. Initial density is a necessary parametric input 
to any solution of the equations of state and is reported, 
for example, by Mader [Ref. 1]. The surface area for the 


reaction zone is somewhat reduced initially because of the 


2 





condensation which occurs in reaction to the departing shock 
wave; subsequent pressure driven expansion of the reaction zone 
will provide a compensating increase in surface area. The 
surface area of the unreacted mass is considered to be a good 
approximation. As the fireball expands, the surface area in- 
creases; the temperature, however, decreases. The rate of 
energy radiation is a linear function of the radiator area, 
but depends on the fourth power of the temperature, as seen 
above in the Stefan-Boltzmann law. Thus, as the fireball ex- 
pands, exiténce decreases. Some argument could be advanced 
that the fireball would encounter additional oxygen from the 
air as it expands, thus helping to maintain the temperature 
through exothermic formation of additional oxides. This has 
been ignored for three reasons: (1) since air is only 20% 
oxygen, for every unit of mass of oxygen taken into the fire- 
ball, four units of mass of neutral (non-reacting) gas are 
also taken in and all the new material taken in must be raised 
to the equilibrium flame temperature, thereby compensating for 
any exothermicity of the oxide formation; (2) except in the 
case of very oxygen deficient or highly metallized explosives, 
formation of additional oxides may not be particularly domi- 
nant because of the decrease in pressure which would accompany 
the expansion; and, (3) the number of oxygen molecules in the 
air immediately surrounding the surface available for reaction 
is small compared to the number of product gas molecules. 
Although afterburn does indeed contribute to overall yield in 


some explosives [Ref. 20], the consideration made is that the 
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temperature at the C-J plane (the detonation temperature) is 
the maximum, and that the exitance computed for that tempera- 
ture and for the surface area of the original explcsive mass 
represents a maximum power level; total maximum energy loss 
through thermal radiation is then established for a time 
period. The numerical value used for this time period varies 
for each explosive, dependent on availiable data. Actual light- 
intensity measurements are available for only some materials; 
the narrow bandwidth of the sensors used to obtain these data 
make them of marginal value. Experimental times of arrival 
of steady overpressure level or other time of significance is 
used where available. The maximum time of significance is in 
all cases determined to be the time to establish the relative 
quiescent peak overpressure which is considered the primary 
damage mechanism in the studies reported by Kinney, Sewell 


and Graham [Ref. 2]. 


D. OPACITY OF THE REACTION ZONE 

Thermal photons emitted interior to the surface do not 
contribute to the radiant flux emanating from the fireball. 
The reason for this is the opacity of the reaction zone. 
The radiation, as photons, is absorbed or scattered before it 
can cross out of the sphere through the reaction zone. Time 
between events for photon-molecule interaction is found using 
the previous relation: 

1 


i= = ; 
S Teo Vv 





v in this case is the velocity of a photon traversing the 


body of gas. AS an approximation, Cor the speed of light in 


2m) 





vacuum is used. o in this case is the cross section of the 
molecules, an admittedly smaller number than the effective 
cross section seen by the radiation as an E-M wave. An order 
of magnitude comparison is again employed; a too-small a, 
which results in a too-large toe is inconsequential to the 
outcome of the comparison since, uSing the average molecular 


10 


diameter of 4x10  m and the n from the thermodynamic equilib- 


rium calculation above, to is found to be on the order of 
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5x10 ““s (4.72x10 ““s in the case of TNT of density 1.60 


gm/em*, from the Mader data). The time for the photon to cross 
the reaction zone of typical 2.0 cm width is ee 00 eS 
Minimum. Thus, it is not probable that a photon from the 
interior would not be scattered or absorbed in the reaction 


zone transit, thereby escaping to contribute to the radiant 


& lec 


E. REFLECTION AND ABSORPTION 

In this model, all thermal radiation emitted by the fire- 
ball is considered to be absorbed by container surface. 
Typical reflectivity profiles for substances used as coating 
material, such as linseed oil based paintS, show a maximum 
reflectivity of about 80%, with values under 50% more common; 
reflectivity is, of course, wavelength dependent. Reference 
21 contains a number of graphs of reflectivity vs. wavelength 
for many coating materials. Radiation emitted from the fire- 
ball will encounter surfaceS tens cf thousands of times per 
millisecond in the container, and be partially reflected and 


partially absorbed at each encounter. After 25 reflections, 
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even wavelengths with a reflectivity as great as .8 would be 
absorbed to less than .5% of their original intensity. About 
.05% of the equatorial cross section of the container is 
composed of the fireball cross section, under the model. 
Absorption in the fireball, with subsequent re~-radiation, can 
be assumed for only some small portion of the radiated energy, 
and is ignored. Atmospheric attenuation within the container 
is likewise ignored because of the small path length involved 


to virtual extinction through surface absorption. 


F. TIME OF RADIATION 

Determination of the actual length of time for which 
thermal radiation losses are a mechanism reducing peak over- 
pressure, and at what rate, involves complexities beyond the 
scope of this investigation. As an absolute maximum, the 
amount of time, called the that this radiant flux reduces the 
internal energy available for overpressure generation is that 
from the time of arrival of the reaction zone at the surface 
of the explosive to the time of establishment of a level cf 
quasi~static overpressure at the surface of the container. 
During this time, the temperature decreases from the C-J tem- 
perature, although by how much is a function of the amount of 
afterburning. 

As a maximum, the total possible exitance is considered 
to be that from a blackbody, radiating at the C-J temperature, 
with surface area equal to that of the unreacted explosive 
mass, for a time, the as described below. Radiant power 


decreases as the fireball grows because the temperature 


rae) 





decreases with volume expansion and exitance is a function of 
the fourth power of the temperature. This power decrease due 
to temperature drop dominates the power increase attributable 
to greater radiator surface area and afterburning where 
present. 

Times of arrival of the initial shock and quasi-static 
overpressure level have been reported for TNT [Ref. 6]. These 
times are each about five milliseconds. If the shock had main- 
tained the C-J detonation velocity for the transit across the 
five meter average distance to the container wall, its time of 
arrival (a little over half a millisecond) would have been 
about one order of magnitude smaller than the measured arrival 
time. Application of the same 10:1 ratio and approximate 
equality of initial shock and quasi-static overpressure arrival 
times to other explosives is an artificiality, although not 
unreasonable for inclusion in the model since the mechanisms 
which cause velocity decay operate on both the initial shock 
and overpressure waves. When no experimental data are avail- 
able, a time of radiation is chosen to be twenty times the 
idealized initial shock transit time as computed from the 
C-J detonation velocity. Where available, actual reported 


arrival times are used. 
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IV. RESULTS 


A. DATA 

Determination of the energy of detonation contained in a 
ten kilogram mass of explosive was rather direct since values 
of g, the specific heat of detonation, are reported in the 
literature, and Q is qxm. The value of gq depends, to some 
extent, on the condition of the explosive; where full data are 
available for more than one set of initial parameters for a 
given explosive, calculations were made for each set. Specific 
volume, v, the reciprocal of density, was determined, since 
the surface area is a constant times the 2/3 power of the spe- 


BD for a constant ten kilogram 


cific volume (A. = 2245 (vo) 
explosive mass. 

As noted above, the explosive energy iS partitioned, with 
some portion going into the initial shock wave. The precise 
partition is highly variable, dependent in part on the density 
of the surrounding medium. Cole [Ref. 22] reports that 533% 
of the explosive energy of TNT is transmitted to the surround- 
ing medium in an underwater explosion. Cook [Ref. 10] examines 
these variances in fluid and condensed media in some detail, 
although no empirical relations are established, per se. For 
purposes of this investigation, an examination of the ratio of 
the total radiant energy emitted to half the total explosive 
energy is deemed appropriate. Division of this (Mxt_) /Q ratio 
by the decimal fraction of Q retained in the explosive after 


departure of the initial shock wave, were it known, would 
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provide the percentage energy loss; it is not considered 
probable that this decimal fraction would be less than .50 
in air. Thus, the maximum radiant energy loss reported in 
the data represents the evaluation of (Mxt_)/.5Q. 

Data and calculation results are contained in Table I. 
Sources of the numerical values used are cited, by reference 
number in parentheses. If no reference is cited, values are 
computed from other entries in the table, using a formula or 
technique described above. The value of q for ALEX is esti- 
mated on the basis of values for similarly aluminized 


explosives. 


B. ANALYSIS 

Limited luminosity data are contained in the literature 
for comparison. Rudlin [Ref. 6, 11 and 19] reported some 
observations on radiant emissions from TNT and pentolite. 
His analysis of data from a series of multi-ton TNT explo- 
Sions, and a scaled comparison to data obtained on a series 
of eight pound sphere detonations are inconclusive due to 
an absence of quantitative measurements of thermal radiative 
power for the smaller charges. In addition, the instrumenta- 
tion used on the smaller masses was band centered toward the 
shorter wavelengths in the visible and near infrared region 
of the spectrum. Rudlin observed that there was substantial 
energy in the infrared not accounted for in his investigation, 
due to the instrumentation red cut-off, at about .7 microns. 
Some data are available for multi-ton TNT events which give 


an integrated, broad bandwidth radiant energy measurement 
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and percentage of explosive yield [Ref. 19]. Unfortunately, 
Pemcmuoce nuclin [Ret. 19], “secaling.....from small charges to 
multi-ton charges, does not, in general, hold." Of interest 
is that luminosity times reported are on the order of a milli- 
second for the smaller charges. Reference 17 aiso contains 
luminosity data. Although no units of power or energy are 
given, luminosity times in a no-afterburn environment (N. 
filled explosion chamber) are also very short, again on the 
order of a millisecond or two for lightly aluminized or non- 
aluminized explosives. 

The times of radiation used in Table I are longer than 
those reported for the small explosive charges and may result 
in an overestimate of the total radiant energy emitted. A 
precise correlation is not possible because of the fact that 
the sensors used to determine the reported experimental times 
were not sensitive to the infrared radiation which makes up 
the majority of the radiant energy emitted. Since maximum 
attenuation effect determination is the object of this study, 
the t. is, within the other model criteria, considered to be 


of an appropriate magnitude. 
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TaBLE I 


Pp, a a. | il 
(gm/om>) |(cal/em)| (u/s) | (°K) 
oxplosive Loss 
Vo Q tL. dik (,>) 
(com>/gm) (Ncal) (ms) (hical/s) 


eed) «| LAGO (7) 2970 (1) 

RDX 08 

PETN 205 
1.54 (7) | 1290 (7)| 6950 (1) | 2940 (1) 

TNT 203 
1.89 (7) | 1480 (7) | 9124 (23)| 2364 (1) 

HMX one 
Vee) | W400 (7) 14784.0.(23)|.2770..(1) 

COMP B Oy 
wee 7) | PASO 47) | 7850 (23)| 2917 (1) 

Tetryl e038 


TATNB or! 
529 eS aes) 0.6041 | 


34 








equwenma ms ee ee a eee. 


joa tol 


Vode i intaam l 
Ceont.) 


is 135 (23) B90 (10))6260 (23/6040 (10) 


044 


DATNB 
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80-20 
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80-20 


45/30/25 
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45/30/25 
RDA-TNT= 
Aluminum 
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V. CONCLUSIONS 


The basic assumption made in previous investigations of 
internal blast [Ref. 2, 3, 4 and 5] that explosion products 
are in thermal equilibrium is shown to be valid. By utiliza- 
tion of a blackbody radiator model, it is shown that explosive 
energy losses through thermal radiation from the surface of 
the explosion fireball are small, ranging from about .5% to 
under .01%. Demonstration of the validity of a second major 
assumption made in these prior investigations, that thermal 
radiative losses could be ignored, is thereby achieved. 

With the recent development of detector systems which 
operate in the infrared, better experimental data on the 
spectral intensity profile of detonating explosives may be 
obtainable to the end of experimentally verifying detonation 
temperatures predicted by thermo-hydrodynamic codes. It is 
probable, however, that the experimental error induced by the 
apparent irreproducibility of explosive results to within a 
tolerance finer than five percent shall remain a major impedi- 
ment to consistently accurate prediction of explosive 


parameters. 
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APPENDIX A 


GLOSSARY OF EXPLOSIVE ACRONYMS 


Aluminized Explosive; 20% aluminum in 802 
Composition B; empirical formula, based on 
100 grams per mole, is 


0 Ik 


Sacee) 2mmee: Wee l3-2.039° 7.333" 
Composition B; mixture of 64% RDX and 36% 
TNT; empirical formula, based on 347.07 
grams per mole, is 


EmGcn Oa7 5. 7 nGSon sO” 


Diamino trinitrobenzene; chemical formula is 


CHEN,O;. 


Cyclotetramethylene tetranitramine; chemical 


formula is C,HON,O,. 


Penta erythritol tetranitrate; chemical 


formula ia C.H_N,0) >; 


Cyclotrimethylene trinitramine; chemical 


Formula is CH O 


BGG 6: 
Triamino trinitrobenzene; chemical formula 


1S CpHEN,O;. 


Trinitrophenyl methyl nitramine; chemical 


formula is CIH.N.O,- 
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TNT 


mriatconal 


Trinitrotoluene; chemical formula is 


CoHL_N20-. 


Aluminized TNT; numbers index percentages 
of TNT and aluminum (e.g., 80-20 is 803 


TNT with 20% aluminum). 
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